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We study a class of models in which N flavors of massless fermions on the half line are 
coupled by an arbitrary orthogonal matrix to iV rotors living on the boundary. Integrating 
out the rotors, we find the exact partition function and Green's functions. We demonstrate 
that the coupling matrix must satisfy a certain rationality constraint, so there is an infinite, 
discrete set of possible coupling matrices. For one particular choice of the coupling matrix, 
this model reproduces the low-energy dynamics of fermions scattering from a magnetic 
monopole. A quick survey of the Green's functions shows that the S-matrix is nonunitary. 
This nonunitarity is present in previous results for the monopole-fermion system, although 
it appears not to have been noted. We indicate how unitarity may be restored by expanding 
the Fock space to include new states that are unavoidably introduced by the boundary 
interaction. 
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1. Introduction 



We consider the action 
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dt 2 + E «*(*)^«( : ^(0,*)Vi(0,t) : - : $(0,^(0,*) : ) 
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i=l i,i=l 



+ [ dxUl(x,t)(d t -d x )i; l (x,t) + i;j(x,t)(d t + d x )ij l (x,t) 



where ^ and ^ are chiral components of Dirac fermion fields, My is an arbitrary real, 
orthogonal matrix, CKj are the rotor degrees of freedom, and the Ii are constants. In 
writing the action, we have split the Dirac fermions into their two single-component chiral 
constituents, ipi and ipi, with the index i labeling flavors. 

The action describes N flavors of Dirac fermions living on the interval < x < L. 
There are no interactions except at x = 0, where the fermion currents are linearly coupled 
to N rotor coordinates. The interaction at x = dynamically couples the left movers ipi 
to the right movers ^ through the rotor coordinates. At x = L, we impose the simple 
reflecting boundary condition ipi(L,t) = e 2 ' K ' lXi tpi(L 1 t) so that the system will be closed 
and self-contained. 

In general, the coupling at x = allows the fermions and the rotors to exchange energy, 
but we will be interested purely in the cases where energy is not exchanged: namely, when 
Ii — > oo or Ii — > 0. If we take all the Ii — > oo, the rotors decouple and a simple reflecting 
boundary condition is imposed on the fermion currents at x = 0. If we take some Ii — > oo 
and the others to zero, then things are much more interesting; we will find that we can 
integrate out the rotators, leaving behind a non-trivial, conformally invariant boundary 
condition on the fermions. The particular boundary condition obtained depends on which 
Ii — > and the value of the matrix My. 

An interesting special case arises when we chose My and Ii so that the action becomes 
essentially equivalent* to one first considered by Polchinski [Q. Polchinski's action captures 
the essential physics of charged fermions scattering from a magnetic monopole in four 
dimensions. 



* Polchinksi uses right moving Weyl fermions on the full line, while we consider right and left 
moving Dirac fermions on the half line. One can convert the half- line Dirac theory into a full-line 
Weyl theory by defining ipi(—x,t) = ipi(x,t) for x > 0, so that left movers at x > get reflected 
into right movers at x < 0. 
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2. The Partition Function 



We use standard bosonization techniques || to rewrite the action in terms of the 
boson fields 



A 



dt 



1 OLiit) 



N 



i=l i,j' = l 

i N r L 



(2.1) 



where $ = 9 4 $ and <&' = The bosons are compact, so $^ = $i + 2n. We fix 

$i(L, £) = 9i, where Q{ = ir(l — 2\), and allow <&j(0,i) to vary with the action. The 
classical equations of motion obtained by varying the action are 

N 

3=1 
N 



&j(0,t) -47r^)d i (t)M ii =0 



$?(M)-$i(M) = 



(2.2) 



Eliminating the rotor coordinates cnj and changing over to the 'rotated' fields <&i(x,t) 



4^1^(0,0 + $ / i (0,*)-C i = 0, 



(2.3) 



where the Cj are constants of integration. 

At this point, let us consider the limit Ii —>■ oo for all i, which should reproduce 
the free fermion theory. The bosons are subject to the boundary conditions $j(0,i) = 
0, t) = 0i and the compactness condition t) = t) + 2-7T, which result in the 
following standard free partition function: 

oo N 

z=n( i -9")" N n 



n=l 



i=l 



(2.4) 



where q = e~^^l L and (3 = inverse temperature. Note the presence of the factor 
Y2^L-oo q^ n+ ^^ , which is simply a sum over winding number. Without the compactness 
condition on this factor would not be present. 
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For the fermion system, we have the boundary conditions ^(0,t) = ^(0, £) and 
ipi(L,t) = e 2n ' lXi tpi(L 1 t). This yields the partition function 

N oo 2 

z =u n (i+? |n - Ai1 ) . (2.5) 

i=l n= — oo 

An application of Jacobi's triple product formula shows that the bosonic and fermionic 
partition functions are indeed equal up to a shift of zero-point energy, verifying our choice 
of boundary conditions for the bosonized action. 

Now let us consider the general case Ii — ■> oo for 1 < i < a and Ii — > for a < i < N . 
We find that 



&!(x,t)-$i(x,t) = 0, 
$i(0,*) = 0, 
^(0,t) = 0, 



JV 



1 < i < N 
l<i<a 
a < i < N 



$i(L, t) = y^ MijOj, 1 < i < N 



rij E Z 



(2.6a) 
(2.66) 
(2.6c) 

(2.6d) 
(2.6e) 



In principle, there should be a constant of integration in equation (|2.6cj ), but we have set 
it to zero because the rotor coordinate will carry an infinite amount of energy otherwise. 

Since equations (|2.6a| ) - (|2.6ej) describe a set of iV uncoupled free bosons, it seems 
that we should be able to compute the partition function as a product of iV independent 
partition functions. The fact that M is an orthogonal matrix means that the canonical 
commutation relations are 



[$i(x, t),$j(y, £)] = [§i{x, t),$j(y, t)} = 4niS i:j S(x - y), 



(2.7) 



so there is no problem with quantizing the system in terms of the The only catch is 
that we must take care in treating the winding modes. When the boundary condition at 
x = is 3^(0, t) = 0, it turns out that no winding mode exists, even though is 
compact. This lead us to the partition function 



71=1 



1 _ q( n ~k) 
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E 9* 
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EU MO 1 



(2- 
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where 

AT 



Z 



M 



|(ni,...n a ) n» = y^^Mjjrij, rij e z|. (2.9) 



The partition function almost looks like that of a free theory, except that the sum over 
winding modes is rather peculiar. Barring the trivial cases a = and a = N, the winding 
mode sum no longer factors into a product of independent winding mode sums. Instead, 
the winding modes of the different boson flavors are coupled together by the matrix M. 

Furthermore, it seems that the winding number sum in the partition function does 
not make sense unless we require M to satisfy a certain rationality property, the details of 
which will be explained in the next section. When M satisfies the rationality property, the 
set Zm is a discrete lattice, so points in Zm are separated by finite gaps. When M does 
not satisfy the rationality property, the set Zm becomes dense on R a . Since we sum over 
the points of Zm in the winding number sum, we find that the partition function diverges 
due to the infinitesimally close spacing of energy levels. 

To flesh out the substance of the preceding remarks, we will work out a couple of 
simple examples: the general N = 2 case and the magnetic monopole for arbitrary N. 



3. N 



Let us consider N = 2, a = 1, 0j = 0, and write M in the form 



M 



1 r 
-r 1 



(3.1) 



This gives us 



Z 



M 



n 



n 



rii + rn2 
Vl + r 2 



; m,n 2 6 Z >. 



(3.2) 



Using some elementary number theory , we see that Zm becomes dense on R unless 
r is rational. For N = 2, the constraint r e Q is the rationality condition for M. Writing 
r = p/l for p,i G Z with gcf(p, I) = 1, we find that Zm ~~ 1 
function is 



=Z, so the partition 



1/ i_^("-^) N 



Em 



(3.3) 



n=l 



m= — oo 
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Using standard techniques [Q, the partition function may be reexpressed in terms of 

4tt 2 

q = e ln <9) ; 

j 2 i 72 00 —1 00 

Z = JP-±Lq-h q is Y[ (l - g n ) (l - <? (2n - 1} ) X] 9 4(p2+ ^ )n2 (3-4) 

n=l n= — 00 

Now we can easily take the limit L — > 00: 

In(Z) - i In(^) + g + O(e-W). (3.5) 

In addition to the standard piece in ln(Z) which scales linearly in the size of the system, 
we see that there is an L-independent term. If we associate the size-independent term with 
the boundary interaction at x = 0, we find that the boundary contributes a temperature- 
independent entropy Sf, = \ ln( p t l ). Equivalently, this means that there are g states 
associated with the boundary, where g = e Sb = y ^ p2 +' 2 , The fact that g need not be an 
integer may seem peculiar, but it is an unavoidable consequence of the way the winding 
modes for different bosons become linked at the boundary. The crucial point is that g > 1, 
giving us a hint that a correct treatment of scattering may need to take into account some 
hidden degree of freedom on the boundary. 



4. The Magnetic Monopole 



Now we take iV > 2, choose a 
I 



N — 1 and 6>, = 0, and set 



M 



1 

V2 
1 

V6 



1 

1 

V6 



2 








y/N(N-l) 
1 

\/jV 



1 



y/N(N-l) 
1 



v /Ar ( Ar " 1 ) 
1 








v /Ar ( Ar " 1 ) 
1 






(N-l) 



1 



/ 

(4.1) 

The last row of M causes the sum of the fermion currents to be coupled to a rotor cxn 
which has / — > 0. The first — 1 rows of M couple to rotors with / — > 00, ensuring 
that differences of fermion currents obey reflecting boundary conditions at x = 0. In 
Polchinski's version [ffj, there is just a single rotor, which correponds to our a^. The 
essence of the model is that we are changing only the boundary condition on the current 
which carries the total U(l) charge. 
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The only bit of work we need to perform is to find a convenient way to make the 
sum over Zm explicit. If we think of the numbers hi as functions hi(rii, . . . , tin) = 
YljLi MijUj, then we see that n^ni, n 2: . . . , tin-i, n,jv + l) = ^{n-y — 1, n 2 — 1, . . . , un-i — 
l,njv) for 1 < i < N — 1, so we can fix un = and just sum over m, n 2i ■ ■ ■ > n^-i- 
Note that we don't care about the value of un since it does not appear in the winding 
mode sum. Furthermore, summing over ni,n2, • • • , n>N-i gives us each possible value of hi 
exactly once; this is true because the matrix obtained by deleting the last row and column 
of M is a non-singular matrix. Finally, we note that J^iLi 1 "I = 12iLi n l ~~ ~k ( J2iL± n i) 2 
Our desired partition function is 

z=n(i- 9 -)" cw ""(i-,c-«)" 1 ± < ,i(E:->?-^(E,~r -.)'). ( 4.2) 

n=l n 4 =— oo 

l<i<JV— 1 

4tt 2 

As in the N = 2 case, we can express Z in terms of q = e ln («) : 

v/2^det(A) n A = H y ; v ; ^ 

l<i< JV-1 

where Ay = — + |<%. det(A) may be evaluated using a standard formula for circulants 
0, yielding det(A) = 2 1 ~ N ±. Taking L — > oo, we find the boundary degeneracy to be 
(7 = Note that # = 1 when N = 2, so we should expect that there are no degrees of 

freedom on the boundary when N = 2. 



5. Green's Functions 

Now we set L = oo and proceed to compute an arbitrary fermionic Green's function 

T: 

T =< Ofe (f . . . $ p (z p )^ h (w x ) . . . ^ (w q )\0 > B , (5.1) 

where we have switched to holomorphic coordinates z = r + ix and imaginary time r = it. 
To distinguish between the presence of tp and ip , we assign Fk = 1 for ^ fc) Ffc = — 1 for 
V>J fe , Gfc = 1 for ^j fe , and Gk = — 1 for ?/4 fc . The notation < ... >b indicates an expectation 
value in the presence of the boundary interaction at x = 0, while < ... > will be used to 
indicate an expectation value for a free theory on the full line with no boundary interaction. 
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As in the finite volume case, we bosonize the system according to the standard corre- 
spondence ipi(z) = e l$ ^ 2 ) and i>i{z) = e - **^^ *. Defining &i(z, z) = Yl!j=i Mij&j{z, z), 
we find that 

d z d s $i(z,z) =0 for3(z)>0 

3>i(z, z) = for 1 < i < a, z = z 
(d z - d B )$i(z, z)=0 for a < i < N, z = z. (5.2) 

We can decompose $i(z, z) into left and right moving fields $f (z) and $f (z), but we 
must remember that the decomposition is not unique since <&i(z,z) = + &f(z) = 

(<lf (z) +C) + ($f (z) -C). Accordingly, let us define constants Coi such that < $f (z) > = 
— < ^f{z) >= Ui. For convenience, define <^>f (z) = $f (z) — (Dj and 4>f{z) = &f{z) + 
Using equation (5.2), we can solve for the left movers in terms of the right movers: 

N 

<Pf(z) = Y J S^(z) 
k=i 

N 

Sik = a j M ji M jk 

i=i 

/ 1, 1 < i < a . . 

ai = \-l, a < i < N ' ^ 

where 4>i = Hj=i M ji4>j- 

The matrix S is manifestly symmetric, and a short calculation shows that it is also 
orthogonal. Furthermore, any matrix which is real, symmetric, and orthogonal may be 
written in the form S%k = — Y2jLi CjMjiMjk with a| = 1 and M real and orthogonal. 
From this, we see that the boundary interaction at x = merely 'rotates' the left movers 
with respect to the right movers, and we may choose the boundary interaction so that any 
particular symmetric and orthogonal matrix S carries out the rotation. 

In bosonized form, equation ( |5.1| ) becomes 

T =< 0\ e - lFl < (Sl) ... e ~ lFp ^ { ~ Zp) e lGl ^i iwi) ...e Gq ^ (Wq) \0 > B . (5.4) 



* To avoid cumbersome anticommuting factors normally appearing in multi-flavor bosonization, 
we simply assume that F is written with the fields in the canonical ordering i\ < 12 < . . -i p , 
ji < 32 < ■ --jq- 
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Using equation (5.3) to write all the &f(z) in terms of the we find that T is simply 

the expectation value of a product of right moving chiral vertex operators. Once V is 
expressed in terms of just right moving fields, we can forget about the presence of the 
boundary interaction and compute as we would in a free theory. Use standard techniques 
for vertex operators H, we get 



N p q 



r^E^/^irp II (*, -„)'-'>'"< II 

a<(3 \<u 

nfi(^-^)- FQGAs — , (5.5) 

where f r = Y^a=i Fa&i a r and g r = X)a=i G\8j xr - Note that f r and (7 r are simply the total 
amount of each flavor present in the ingoing and outgoing states, respectively. 
For the magnetic monopole, SV, = — Sij + giving us 

N p q 

r m ono P ole = e iKU, l[8 fr+griK f[ (z a -Zp) FaFp5 ^ j[ {w X - w v ) GxGv&i »" 

r=l ac, 0=1 X,t/=1 

a<f3 \<u 

nn(4-^) FA(s,airi) , (5.6) 

a=l A=l 

where k = -j| Ylr=i 9 r an< ^ 00 = ^2r=i °° r - ^ ur resu fts for the monopole Green's functions 
agree with the original calculation by Polchinski , in which he evaluated the path integral 
by integrating out the fermions. The parameter 10 in our result corresponds to the vacuum 
instanton angle 6 in Polchinski's result. The path integral calculation gives results which 
are integrated over all values of 9, from which Polchinski extracted the ^-vacuum Green's 
functions by cluster decomposition arguments. In our bosonized solution, we automatically 
obtain the Green's functions in a fixed ^-vacuum sector. 



6. A Unitarity Paradox 

Now we would like to point out that the Green's functions of equation (5.6) produce 
an apparently nonunitary S-matrix for N > 2*. Let us consider Green's functions where 

* The unitarity problem is not particular to our treatment of the monopole-fermion system; 
the results obtained in Q have exactly the same problem, although it seems not to have been 
commented on until now. 
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we have some fixed combination of left moving operators and any number of right moving 
operators. The idea here is that we fix an ingoing state created by left moving operators, 
allow scattering to take place, and then look at the overlap with outgoing states made of 
arbitrary combinations of right movers. Since we have fixed the left movers, the integer 
parameters G\ and g r are fixed, while F a and f r are free to vary. Since k = X)r=i 9ri 
we see that k is fixed, but not necessarily integral for N > 2. Looking at equation (5.6), 
we see that the Green's function will vanish unless f r + g r = k for all r, leading us to 
the following conclusions: for N = 2, every choice of ingoing state has overlap with some 
outgoing states, and there is no problem with unitarity. In fact, the interaction at N = 2 
simply swaps the two flavors upon reflection from the boundary. For N > 2, we find that 
any ingoing state with a non-integer k (i.e., a state whose total charge is not a multiple of 
N/2) has zero overlap with all possible outgoing states, violating probability conservation 
and unitarity. 

As a first step toward resolving the unitarity problem, let us think about the S-matrix 
in terms of bosons. The advantage of the bosonized form is that we can think of the 
scattering directly in terms of operators: by virtue of equation (5.3), we can think of the 
boundary S matrix as an operator which maps the chiral, left moving 'in' state operators 
into the chiral, right moving 'out' state operators. Explicitly, we have S4>^{z)S~ l = 
Y^j=\ Sij(pf(z), where S is the Hilbert space S-matrix operator and we regard (pf(z) and 
4>f"(z) as operators in a free theory without a boundary interaction. 

Formally, it is easy to show that S acts on the bosonic Fock space in a unitary way. 
Let | in > and \in' > denote two arbitrary ingoing states formed by acting on the vacuum 
|0 > with arbitrary combinations of operators (j)f(z). To demonstrate unitarity, we need 
to show that < in'\S'S\in >=< in'\in >. If we let S and act on the left moving fields 
in | in > and < in'|, the left movers turn into right movers multiplied by factors of Sij. 
Since we are now dealing with free field theory, the inner product may be evaluated in 
terms of two-point functions by Wick's theorem. The orthogonality of Sij implies that 
<0|#(z)#£#(u;)|0> = Em=i^^<0|^(z)^H|0> = <O|0f(^H|O>, 
which finishes the proof that S is unitary. 

Now let us consider what happens when S acts on fermionic states. Setting N = A 
for simplicity, we find that 

sm^s- 1 = se^^s- 1 = e -i*f(-) e i* 2 R (-) e i*3 R w e i<w 

s^s- 1 = se^^s- 1 = e i*f w e -i*?w e i*?w e i*f ( 6 .i) 
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with similar expressions for the other flavors. Every ingoing fermion operator scatters 
into a product of vertex operators carrying half-integer charges of each flavor. Since the 
expectation value of a product of vertex operators vanishes unless there is a net charge of 
zero for each flavor, we see that ingoing states with an odd number of fermions have no 
overlap with any outgoing states made of fermions. This suggests a simple solution to our 
unitarity problem: we should expand our Fock space and allow outgoing states to contain 
vertex operators carrying half-integer charges. The new operators are in fact nothing but 
spin operators for our Dirac fermions. The complete Fock space will be generated from 
fermion operators and spin operators, and the boundary interaction will couple the two 
kinds of operators together. We can also think of the Fock space as being built from fermion 
operators acting on two kinds of vacuum states: one is the ordinary Raymond vacuum, and 
the other is the Neveu-Schwarz vacuum created by a spin operator acting on the Raymond 
vacuum. When an even number of fermions scatter from the boundary, we remain in the 
same vacuum. When an odd number of fermions scatter from the boundary, the outgoing 
states emerge in the opposite vacuum from the ingoing states. In some sense, we can think 
of the boundary as having a two-fold degeneracy corresponding to the two possible vacua. 
We plan to present a more detailed exposition of these ideas in a forthcoming work. 

7. Conclusions 

In this letter we have generalized the monopole-fermion system to a class of conformal 
systems with conformal boundary interaction parameterized by an orthogonal matrix. 
Abelian bosonization allows completely explicit computation of the partition function and 
Green's functions, which agree with previous results obtained for the monopole-fermion 
system. We point out that the scattering appears to be nonunitary when N > 2, and 
suggest how this may be corrected by recognizing the presence of additional states of the 
boundary. The partition function tells us that the boundary degeneracy is g = \J~^, so 
extra states for the boundary are required precisely when g > 1. 
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